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Direction Finding 

This invention relates to direction finding, and more particularly to direction finding using 
radio techniques. 

Direction finding by radio is known. In a direction finding system having a number of 
5 separate antennas, radio wavefronts reach the antennas with different delays. Assuming 
the wavefronts have narrow bandwidth, these delays give rise solely to relative phase shifts. 
For determining the directions or bearings of M-1 radio frequency (RF) emitters (M is a 
positive integer), a typical direction finding system employs M antennas, a respective 
receiver for each antenna, and processing circuitry for each receiver. The processing 

10 circuitry implements a discrete Fourier transform to divide its respective signal into 
frequency bins or channels. Signals are then combined in pairs, and for each channel a 
respective covariance matrix is constructed from which emitter bearings are estimated. This 
approach suffers from the problem of requiring as many receivers as there are antennas, 
receivers being expensive and bulky. It is unsuitable for example for man-portable 

15 equipment or for mounting in small aircraft. 

To get round the multiple receiver problem, it has been proposed to use a single port 
receiver, spatial spectrum estimation technology, and a weight perturbation algorithm to 
obtain the covariance matrix. See Zhao Yimin, "A Single Port Receiver Spatial Spectrum 
Estimate DF System", 0-7803-3216-4/96 IEEE 1996. This is however a relatively complex 
20 approach to the problem. 

It has also been proposed to combine RF antenna signals in beamformers and use a single 
receiver connected to successive beamformers via a multipole switch, one pole per 
beamformer. See CMS See, "High Resolution DF with a Single Channel Receiver", 0- 
7803-7011-2/01 IEEE 2001. A single spatial covariance matrix is formed, and the bearings 
25 or directions of up to M-1 emitters can be estimated, where M is the number of antennas. 
This requires M 2 beamformers, another source of expense and bulk. 

Another solution is adopted in the Rohde & Shwartz DDF 195 instrument, which combines 
pairs of RF antenna signals with each of four relative phase shifts inserted between them in 
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succession by means of switches This requires only a single receiver, and uses one 
antenna as a reference antenna, combining its output with that of other antennas in turn 
with multiple switched phase shifts. However, the method is for estimating the bearing of a 
single emitter only. A patents search has indicated that the following patent documents 
5 EP4551 02, DE401 4407, DE3636630, DE1 9529271 and DE2723746 are related to direction 
finding. 

US patent application no. US 2002/0190902 A1 describes sampling RF signals from M 
antennas at a sampling rate equal to or greater than the signal bandwidth multiplied by 2M. 
This is followed by Fourier transformation of resulting signal samples to provide spectra. 
10 Direction finding is then based on line configuration in the spectra and associated phase 
and amplitude data. It employs a multiplicity of directional antennas to provide 
beamforming. 

US patent no. 4,486,757 to Ghose et al. discloses a direction finder having two separated 
antennas for reception of respective signals from a remote source, phase shifting one signal 
15 by 90 degrees relative to the other and nulling one signal using an error correction loop, 
the operation of the loop is monitored and used to calculate a bearing for the remote 
source. It does not appear to use receivers. 

Direction finders of various kinds are also disclosed in US patent no. 4,489,327 to Eastwell, 
20 British patent no. 1 ,576,61 6 and US patent application no. 2002/0008656 to Landt. 

It is an object of the present invention to provide an alternative form of direction finding. 

The present invention provides a direction finding system incorporating a plurality of 
antennas characterised in that the system also includes: 

25 a) means for determining individual antenna signal strengths; 

b) combining means for determining combined antenna signal strengths by forming 
combinations of first and second antenna signals derived from different antennas, 
wherein the second antenna signals are in two sets with signals in one set having a 
non-zero phase difference relative to signals the other set; and 
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c) means for determining at least one emitter bearing from antenna signal strengths 

The invention provides the advantage of requiring only a single receiver when successive 
signal strengths are determined in successive steps, which reduces cost and bulk and still 
provides a viable direction finding technique. 

5 The means for determining emitter bearing may be arranged to derive covariance matrix 
elements from antenna signal strengths and to determine emitter bearing therefrom. It may 
alternatively be arranged to derive a relationship between antenna signal strengths and 
emitter bearing and to determine emitter bearing therefrom. 

The relative phase difference may be in the range 30 to 120 degrees, preferably 
10 substantially 90 degrees. 

The combining means may be arranged to combine antenna signals with equal gain 
magnitude and with equal and unequal phase. It may incorporate phase shifting means 
switchable into and out of an antenna signal path, and an adder having two inputs both 
switchably connected to individual signal paths extending to respective antennas. 

15 The means for determining individual antenna signal strengths may comprise a first 
multipole switch having input poles connected to receive signals from respective antennas; 
the combining means may incorporate a second multipole switch having input poles 
connected to receive signals from respective antennas and a third multipole switch for 
switching phase shifting means into and out of an antenna signal path extending via the 

20 second multipole switch; and the combining means may also incorporate adding means for 
combining signals, the adding means being arranged to add an antenna signal in a first 
signal path extending via the first multipole switch to another antenna signal in a second 
signal path extending via the second and third multipole switches. 

In another aspect, the present invention provides a method of direction finding using a 
25 plurality of antennas characterised in that the method incorporates: 

a) determining individual antenna signal strengths; 
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b) determining combined antenna signal strengths by forming combinations of first and 
second antenna signals derived from different antennas, wherein the second antenna 
signals are in two sets with signals in one set having a non-zero phase difference 
relative to signals the other set; and 

5 c) means for determining at least one emitter bearing from antenna signal strengths. 

Emitter bearing may be determined by deriving covariance matrix elements from antenna 
signal strengths and determining emitter bearing therefrom. It may alternatively be 
determined by deriving a relationship between antenna signal strengths and emitter bearing 
and determining emitter bearing therefrom. 

10 The relative phase difference may be in the range 30 to 120 degrees, preferably 
substantially 90 degrees, and successive signal strengths may be determined in successive 
steps. 

The step of forming combined antenna signal strengths combines antenna signals with 
equal gain magnitude and with equal and unequal phase. It may include switching phase 
15 shifting means into and out of an antenna signal path, and adding signals in signal paths 
extending switchably to different antennas. 

The step of determining individual antenna signal strengths may comprise switching signals 
from antennas via a first path. The step of forming combined antenna signal, strengths may 
incorporate: 

20 a) switching signals from antennas via a first path for combining; 

b) switching signals from antennas via a switch selectable second path or a third path for 
combining, the third path being arranged to phase shift antenna signals therein relative 
to antenna signals in the second path; and 

c) adding a first path antenna signal to second and third path antenna signals individually. 

25 

In order that the invention might be more fully understood, embodiments of the invention will 
now be described, by way of example only, with reference to the accompanying drawings, 
in which: 
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is a schematic block diagram of a direction finding system of the invention; 

is a graph illustrating direction finding results obtained using the invention; 

is a graph illustrating direction finding results obtained using a single receiver 
prior art system comparable with the invention; and 

is a graph illustrating direction finding results obtained using a multiple receiver 
prior art system; 

Referring to Figure 1, a direction finding system of the invention is indicated generally by 
10. As illustrated for the purposes of this example, the system 10 incorporates four 
antennas 12 each of which is conventional, i.e. omnidirectional. In general, as many 

10 antennas may be used as are required to detect a desired number of emitters, i.e. M 
antennas for M-1 emitters. Signals pass from the. antennas 12 via respective buffer 
amplifiers 14 to first and second multipole switches SW1 and SW2, the amplifiers 12 being 
connected to respective input poles b, c, d and e of both switches. The multipole switches 
SW1 and SW2 have respective movable contacts C1 , C2 which allow any of the associated 

15 input poles b to d in each case to be connected to respective output poles f. 



Figure 1 
Figure 2 
Figure 3 

5 Figure 4 



The output pole f of the first switch SW1 is connected to a first phase shifter P1 introducing 
a phase shift of fa t and the output pole f of the second switch SW2 is connected to second 
and third phase shifters P2 and P3 introducing phase shifts of <f> 2 and (fo respectively. In the 
present example, where an ideal situation is envisaged, all three phase shifters P1 to P3 

20 have a gain of unity: the first and second phase shifts fa and <j> 2 are zero (and the 
associated phase shifters P1 and P2 could be removed and replaced by connections), and 
the third phase shift fa is 90 degrees. However, non-ideal gains/phases can be 
accommodated by calibration (as will be described later). Satisfactory performance has 
been demonstrated with errors up to plus or minus 60 degrees: in other words the phase 

25 difference <J> 3 introduced by the third phase shifter P3 may be anything in the range ±30 to 
±120 degrees, with simultaneous gain discrepancy of up to 3 dB compared to unity. Larger 
discrepancies result in progressive deterioration of .performance. This implies the third 
phase shifter P3 providing for <j) 3 to. be preferably in the range 30 to 120 degrees different to 
<|>2 but fa is unrestricted. It is also acceptable for phase shift and gain to vary with frequency, 

30 provided the phase shift is reasonably close to 90 degrees and phase shifter gains are 
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reasonably close to 1.0, and that these parameters are stable and predetermined functions 
of frequency. 

An output signal from the first phase shifter P1 passes to a first input 16a of an. adder 18. 
Output signals from the second and third phase shifters P2 and P3 pass to first and second 
5 input poles p and q respectively of a third multipole switch SW3, which has a third input pole 
r connected to ground, The third switch SW3 has a movable contact C3 which allows any of 
the associated input poles p, q and r to be connected to an output pole z, which in turn is 
connected to a second input 1 6b of the adder 1 8. 

Figure 1 is a somewhat idealised drawing in which all components are assumed to be 
10 perfect and (where appropriate) matched. In practice, in a non-ideal situation, components 
may need to be trimmed or impedance matched to counteract unwanted effects by inserting 
additional circuitry, this is well known in the art of electronics and will not be described. 

The adder 18 has an output 20 connected to a series-connected chain of elements 
consisting of a receiver 22, an analogue to digital converter (ADC) 24, a discrete Fourier 

15 transformer 26 and % a digital signal processor (DSP) 28. The digital elements 24 to 28 may 
be replaced by equivalent analogue processing if desired. The transformer 26 provides a 
discrete Fourier transform (DFT) of length N points with windowing. This is a well-known 
calculation and will not be described: see e.g. A V Oppenheim and R W Schafer with J.R. 
Buck, "Discrete Time Signal Processing", Prentice Hall, Englewood Cliffs, NJ. USA, 1999. 

20 The DFT is advantageously computed by fast Fourier transform (FFT), but this is not ' 
essential. The choice N=1 corresponds to the trivial case equivalent to no discrete Fourier 
transform present. 

The direction finding system 10 operates as follows. By appropriate choice of the positions 
of the movable switch contacts C1 to C3, any antenna signal may be added to any other 
25 antenna signal, or selected alone by adding it to a zero signal at grounded input r of third 
switch SW3. Assuming that the three phase shifters P1 to P3 have equal gain magnitudes, 
antenna signals at the adder inputs 16a and 16b have a relative phase shift between them 
equal to <J> r <t>2 or 4h-fa according to whether the third switch movable contact C3 is 
connected to its first or second input pole p or q respectively. In the ideal case, = = 0 



WO 2005/010550 



PCT/GB2004/002902 



7 

and <t> 3 = 90 degrees. In this case, the relative phase shift is 0 or 90 degrees (ignoring sign) 
between adder input signals according to whether the third switch movable contact C3 is 
connected to its first or second input pole p or q. 

The receiver 22 has a front-end bandpass filter (not shown) to attenuate unwanted signals 
outside a frequency band in which the antenna signals appear. It processes the adder 
output signal and converts it to a complex output in a base band convenient for sampling 
and analogue to digital conversion at 24. The sampling rate of the ADC 24 must be 
sufficiently high to avoid aliasing, i.e. it must exceed the bandwidth of the receiver front-end 
bandpass filter. If, for example, the front-end bandwidth is 25 MHz, and a 512-point DFT is 
implemented by the transformer 26, the ADC sampling rate would be at least 27.5MHz, and 
it would take approximately 20 ^s to acquire 512 samples for the transformer 26. The digital 
signal from the ADC 24 is converted by the transformer 26 into a spectrum in terms of a set 
of frequencies expressed as bins or channels of finite width each with an associated ' 
magnitude. The frequency channels would be approximately 100 kHz wide for. a 512 point 
DFT with' sampling at 27.5MHz. A DFT is implemented with windowing to reduce leakage 
between channels. The options for choice of window include a rectangular window 
equivalent to no window. 

Optionally, selection of settings of movable contacts C1 to C3 may be made in such a way 
that not all antennas contribute. This gives faster processing, but possibly less accuracy, 
and it reduces the maximum number of emitters that can be detected. 

As has been mentioned, the various possible settings of the switches SW1, SW2 and SW3 
allow the receiver 22 to input either the signal from any individual antenna 12, or a sum of 
relatively phase shifted signals from any pair of antennas 12. The switches SW1 to SW3 
may be operated to give a random or pseudo-random selection of antenna signals to avoid 
possible deleterious effects with particular signals. The output of the receiver 22 is sampled 
by the ADC 24 and processed by the transformer 26, which computes the N-point 
windowed DFT of a block of N consecutive samples from the ADC. This is a well known 
procedure and will not be described in detail. The output of the transformer 26 comprises N 
frequency domain samples, i.e. magnitudes of the contents of the N frequency bins. For the 

general or nth frequency bin, the frequency domain sample is denoted by S n , where n is a 
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frequency domain index in the range 0 to N-1. S£ is defined as the discrete Fourier 

transformer output with frequency index n if the kth antenna is connected through the first 
switch SW1, the third switch SW3 is connected to OV at r, and the gain of the entire path 
from the kth antenna to the receiver input is unity. 

5 The process of direction finding consists of estimating the bearing or angle of incidence or 
angle of arrival of one or more signals received by the antennas 12. It is carried out for one 

or more frequency bins by the DSP 28, which processes frequency domain samples S n . ■ 

A complex gain constant Gi k is now defined to represent both gain and phase shift applied 
to a radio signal in a path through the kth antenna 12 (k = 1 to M), the first switch SW1 , the 
10 first phase shifter P1 and the adder 20 to the receiver 22. Similarly, complex gain constants 
GaAk and G 2 Bk are now defined to represent both gain and phase shift applied to a radio 
signal in paths to the receiver 22 through the kth antenna 12 via {inter alia) the second and 
third phase shifters P2 and P3 respectively. 

With the first switch SW1 connected to the buffer amplifier 14 of the kth antenna and the 
15 third switch SW3 connected to 0V at r, the transformer output with frequency index n is the 

product G 1k S£ . This transformer output has a mean squared value or power P kk associated 

with the kth antenna and given by: 

P k K=E{|G 1k Sj? | 2 } = |G 1k | 2 E{|5n 2 } (1) 

where E{..} is the mathematical expectation operator and |.,.| (as in e.g. |S£ |) represents a 

20 modulus. The DSP 28 computes a measurement (or "estimate") of P kk which either equals 

a single value of S n or equals an average or weighted average of several values of 5 n , 
these values being obtained from the discrete Fourier transforms of respective blocks of 
data collected with switch contacts C1 to C3 set to the appropriate positions. These blocks 
of data may be overlapping or non-overlapping. For clarity in the following explanation, the 
25 notation Pw< will denote the measurement of P kk obtained as described above. 
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The first switch SW1 is connected to the buffer amplifier 14 of each antenna 12 in turn, i.e. 
the antenna index k goes from 1 to M where M is the number of antennas, and the third 
switch SW3 remains connected to 0V at r. The transformer output power P kk is measured in 



5 The first switch SW1 is now connected to the kth antenna, the second switch SW2 is 
connected to the mth antenna, and the second switch signal path with gain G 2Am is selected 
by connecting the third switch movable contact C3 to its first input . p. The power P kmA 
associated with gain G^m at the transformer output with frequency bin index n is then 
measured, and it is given by: 



where Re{..} represents "real part of" and the asterisk "*" a complex conjugate. P kmA is 
measured for all possible pairings of different antennas 12, of which there are M C 2 pairs. 
The value of P kmA is measured by the DSP 28 in the same way as it measures P kk , as 
15 previously described, and the same convention is adopted that the notation P kmA is used in 
what follows to refer to the measurement. 

In the same way, the procedure associated with Equations (2) and (3) is now repeated, 
except that the second switch signal path with gain G 2Bm is now selected by connecting the 
third switch movable contact C3 to its second input q to implement phase shift instead of 
20 fa. The power P kmB associated with gain G 2Bm at the transformer output with frequency index 
n is then measured for each antenna pairing, and following equivalents of Equations (2) and 
(3) may be generated by replacing index A with index B. 



each case. 



10 P kmA =E{|G 1k S£ +G 2Am S£ | 2 } 



(2) 



i.e. P kmA = |G lk | 2 E{| S£ | 2 } + KWE{| S£ | 2 } + 2Re{ G 1k G* 2Am E{ S£S£ * }} 



(3) 



P kmB = E{|G, k S£ + G 2Bm S£| 2 } 



(4) 



i.e. P kmB = |G 1k | 2 E{| S£ | 2 } + |G 2Bm | 2 E{| S£ | 2 } + 2Re{ G 1k G* 2Bm E{ S£S£ * }} 



(5) 
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In order to carry out direction finding, the gains G 1k , G^m and G 2B m must first be determined 
at every frequency of interest, i.e. at the centre frequencies of those frequency bins defined 
by the DFT operation in the transformer. 26 which are associated with emitters. In one 
approach, this may be done by making electrical measurements of the complex gain of the 
5 individual components of the system, or of groups of components. Such a procedure is well 
known in the art of electronics and will not be described. 

In an alternative embodiment, the determination of the gains may advantageously be 
carried out by a procedure as follows. A radio emitter is set up in a known location, and the 
direction finding system 10 is then operated to make the measurements described in 

10 Equations (1) to (5), for each combination of pairs of antennas, by setting switch contacts 
C1 to C3 to the appropriate positions. Using the known values of the emitter power, 
bearing and distance from the antennas, the quantities E{| S£ | 2 }, E{| | 2 }, and E{ S£S? n * } 
in Equations (1) to (5) may be calculated at each frequency of interest. The only unknown 
quantities in Equations (1) to (5) are then the values of the gains G lkf GaAm and G 2B m, which 

15 may therefore be determined using a procedure known as non-linear optimisation (see e.g. 
M S Bazaraa, H D Sherali, C M Shetty, "Nonlinear Programming - Theory and Algorithms", 
Wiley, New York, 2 nd Ed, 1993). This is well known in the art of scientific computing and will 
not be described. • The accuracy of the estimates of the gains G 1k , G 2A m and G 2B m will 
usually be improved by moving the emitter to one or more further known positions, 

20 repeating the collection of measurements described in Equations (1) to (5), and using the 
entire set of measurements in the non-linear optimisation procedure. 

Once the gains have been determined, one procedure for carrying out direction finding is as 
follows. The value of P kk is measured for each antenna 12 in turn and the value of E {\S£ | 2 } 
is then estimated using Equation (1) and the previously determined value of G 1k as: 

25 E{|5^| 2 } = PKK/|G lk | 2 ' (6) 

Similarly the values of P kmA are measured for each combination of antennas 12 in turn. 
Using the previously determined values of G 1kl G^m, together with the values of E{|S£ | 2 } ' 
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and E{|S£ | 2 } computed above, and Equation (3), a quantity x is now computed for each 
antenna pairing from: 

x= (1/2)(P kmA - |G lk | 2 E{|5n 2 }--|G 2Am | 2 E{|5^| 2 }) (7) 
i.e. x = Re{ G 1k G* 2Am E{S£S£ *}} (8) 

5 In the same way, the values of P krnB are measured for each pair of antennas 12 in turn. 

Using the previously determined values of Gi k , G 2Bmi the values of E{|S£ | 2 } and E{|S£ | 2 } 

computed above, and Equation (5), a quantity y is now computed for each antenna pairing 
from: 

y= (i/sxp^bHG^ ( 9 ) 

10 i.e. y = Re{G 1k G* 2 BmE{ SJX * }} (1 0 ) 

The known complex value G 1k G* 2Am is now written as c + jd, and that of G 1k G* 2Bm as e + jf, 
where j is the square root of -1. The next step is to compute the unknown complex value 

E i s k s m *} written as a + jb. Rewriting Equations (8) and (10) in terms of a to e and j: 

x = Re{(c+jd)(a+jb)} = ca - db (11) 

15 and y = R e {(e+jf)(a+jb)} = ea - fb (1 2) 

Equations (11) and (12) are two simultaneous equations in two unknowns which are solved 
by standard methods to give the required values a and b, which in turn give E{S£S£ *}. 

As an example, for G 1k = 1, G 2A m = 1, and G 2Bm = j for $ = 90 degrees, then c = 1, d = 0, 
e = 0 and f = 1 , and a = x, b = -y. 
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The procedure associated with Equations (6) to (12) is repeated for each chosen pair of 
antenna index values k and m, and the corresponding value of E{S£s£ *} is computed in 
each case. These two procedures yield a set of values E{|S£) 2 } for k = 1 to M and 
*} for k = 1 to M-1 and m = k+1 to M which are known collectively as "covariance 
5 terms". When arranged in a square array with E{|S£ | 2 } terms on the array diagonal (row k 

and column k) and each E{ * } term at a respective row position k and column position 
m, this set of values is known as the "spatial covariance matrix". The entire process may be 
carried out independently for each DFT frequency bin, i.e. each value of transformer output 
index n. 

10 Using some or all of the covariance terms in the spatial covariance matrix, the bearing (or 
angle of incidence or arrival) of one or more received signals may be estimated using 
standard techniques: see e.g. H L Van Trees "Optimum Array Processing" (part IV of 
"Detection, Estimation and Modulation Theory"), Wiley, New York, 2002, which discloses for 
example the MUSIC algorithm and least squares fitting. This is well known in the art of 

15 direction finding and will not be described. 

The foregoing example invention has been described in terms of measuring antenna signal 
powers. The invention may however be implemented with any power-related measurement 
of antenna signal strength: in this connection antenna signal strength is defined as a 
measurement of the discrete Fourier transformer output with frequency index n in the above 
20 example, when the kth antenna is selected by switch SW1 and the input pole of SW3 is 
connected to ground, from which the mean squared value or power P kk used in Equation (1) 
can be derived. Examples of antenna signal strength include mean squared voltage or 
current, or root mean squared voltage or current, in addition to signal power itself. 

25 Similarly, combined antenna signal strength is defined as the discrete Fourier transformer 
output with frequency index n, when the receiver input is the sum of the signal from the kth 
antenna connected through switch SW1 and the signal from the mth antenna connected 
through switch SW2, phase shifter P2 or P3, and switch SW3, from which the mean 
squared value or power P km used in Equations (2),(3) can be derived. 
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It is possible to compute the bearing of a single emitter at the frequency corresponding to 
index k in Equations (1) to (12), without computing the spatial covariance matrix, by an 
alternative procedure which will now be described. It has been found that this procedure 
can be faster than that which computes the spatial covariance matrix because it requires 
5 fewer measurements, but it is less accurate. 

Two antennas are selected, numbered k and m, and the quantities P kmA and P km8 then are 
measured as described earlier. Using standard physical analysis (see e.g. H L Van Trees 
"Optimum Array Processing" (part IV of "Detection, Estimation and Modulation Theory"), 
Wiley, New York, 2002) for the case of a single plane-wave received radio signal, it may be 
10 shown that the values of E{|S£) 2 } and E{|S£| 2 } are equal: for convenience this value is 
denoted as V, i.e. E{| S* | 2 } = E{| S% | 2 } = V. Similarly, it may be shown that E{ S£s£ * } = V 
exp (j0), where 6 is a mathematical function of the difference between the bearing of the 
incoming plane wave and the bearing of the line joining the antennas k and m. Equation (3) 
may therefore be rearranged to give: 

15 P kmA = V(|G 1k | 2 + |G 2Am | 2 + 2Re{ G 1k G* 2Am exp (j9)} ) (13) 
and Equation (5) may similarly be rearranged to give: 

Pkms = V(|G 1k | 2 + |G 2Bm | 2 + 2Re{ G 1k G* 2Bm exp (j9)} ) (1 4) 

Next, antenna k and a further antenna, numbered p, are selected. The quantity P kpA is then 
defined as equivalent to P^a but for antennas k and p. P kpA is measured, and, by a similar 
20 analysis to that for Equations (13) and (14), may be shown to. give: 

P kP A = V(|G lk [ 2 + |G2Ap| 2 + 2Re{ G 1k G* 2Ap exp (\\ V )) ) (1 5 ) 

where y is a mathematical function of the difference between the bearing of the- incoming 
plane wave and the bearing of the line joining the antennas k and p. The three Equations 
(13) to (15) contain three unknown quantities, 9, y, and V. 6 and y are related to each 
25 other because the angle between the line joining antennas k and m and the line joining 
antennas k and p is known. Given the three measurements P kmA , P kmB and P kpA , the 
unknown bearing of the incoming plane wave may be determined, for example by using 
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non-linear optimisation, as previously mentioned. It is not necessary to estimate the value 
of V. 



5 



10 



15 



20 



The principles set out above may be used in the case of more than one plane wave arriving 
from different bearings, to formulate equations to replace Equations (13) to (15). Additional 
equations are formulated in the same way for one or more further pairs of antennas (for 
example antenna p and antenna q), and the unknown bearings may be determined, again 
by using non-linear optimisation. 

Signal power received by antennas 12 may fluctuate: because of this it is advantageous to 
carry out averaging in connection with the measured quantities P kk , P kmA and P k(nB defined 
above. It is convenient to define a "commutation cycle" as consisting of a cycle of collecting 
data as described above from each of the settings b to e, p to r, of the three switches SW1 , 
SW2 and SW3 required to determine P kk , P kmA and P kmB for one value of the frequency bin 
index n. Further data is then collected using additional commutation cycles, each cycle 
giving a respective set of values of P kk , p^ and P kmB . The resulting multiple values of P kk , 
P kmA and P kmB are then used to provide average values of each. These averages are then 
used in Equations (1) to (10) above. 



The rate at which commutation cycles are carried out is referred to herein as the 
"commutation cycle rate". One or more of the signals received by the antennas 12 may 
contain periodic fluctuations, if for example data symbols are. carried by a signal at a 
particular rate. If the period of the fluctuations is equal of close to a multiple or sub-multiple 
of the commutation cycle rate, then covariance terms derived as described above may be 
subject to consistent or systematic errors (known as "biases"). These may result in errors in 
estimated angles of incidence. To avoid this problem, the order in which the commutation 
(switch setting) is performed within each commutation cycle may advantageously be varied 
25 in a random or pseudo-random sequence between successive cycles. 

Equations (1) to (15) above may be re-evaluated for other values of the frequency bin index 
n, e.g. adjacent frequency bands corresponding to index n-1 and n+1. Output bearing 
estimates are then obtained for various values of n and are combined to give an estimate 
for which error limits can be calculated. 
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A simulation was made of a direction finding system 10 of the invention having four 
antennas 12 in a square array and using the approach of calculating the spatial covariance 
matrix. The simulation envisaged two emitters with identical carrier frequencies to be 
located in bearing in the plane of the antenna array. The emitters were treated as both 
transmitting at 6.25ksymbol/sec with a modulation type of QPSK, as described by R E 
Ziemer and R L Peterson, "Introduction to Digital Communication", Maxwell Macmillan 
International, New York, 1992. They were at respective bearings or angles of arrival of 30° 
and 70° at the antenna array relative to a predefined reference direction. The signal to noise 
ratio (emitted signal power divided by total noise power over 25.6MHz bandwidth) was 
assumed to be high. The sampling rate was 25.6MHz and a block of 4096 samples were 
collected, a Hamming window was applied (see Oppenheim et al. mentioned above), and 
the result was processed by discrete Fourier transform. The sample collection time was 
0.16 milliseconds, which corresponds to one emitter symbol period. If only one block of 
samples were to be processed, results are poor, because signals are highly correlated over 
the sample collection time. First and second blocks of samples were collected with an 
intervening time interval corresponding to at least three symbols to avoid this. Emitters of 
interest often contain a root raised cosine (RRC) pulse shaping filter which is shorter than 
this intervening time interval, in which case there is no correlation between signals in the 
sample blocks. The overall covariance matrix was then derived using parameters obtained 
using Equations (1) to (10) for both data blocks and averaging. Matrix elements were 
processed using the known MUSIC algorithm (see Van Trees mentioned above) to 
determine emitter bearings. This algorithm is neither the only nor even the best algorithm 
for this purpose, but it is simple and often suggested for that reason. 

Figure 2 is a graph of MUSIC spectrum against angle of arrival (AoA) derived from 
covariance matrix elements obtained using the method of the invention. It shows 
reasonably well-defined peaks 40 and 42 at angles of arrival of 30° and 70°, showing that 
the simulated emitters have been located. 

For comparison, a trial was made of the single receiver technique disclosed by C M S See, 
"High Resolution DF with a Single Channel Receiver", 0-7803-7011-2/01 IEEE 2001. Using 
the same parameters and MUSIC algorithm processing, this technique gave the results 
shown in Figure 3. Here the emitter at 30 degrees has not been resolved, and the 70 
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degree emitter has been resolved at 50, but with a peak height ~15x smaller than the peak 
42 (NB Figures 2 and 3 have axes with different scales). 

For completeness the conventional method (which requires one receiver per antenna) was 
also simulated using the same parameters and MUSIC algorithm processing. This method 
5 gave the results shown in Figure 4. It shows very well defined peaks 60 and 62 at angles of 
arrival of 30° and 70°, showing that the simulated emitters have been located with better 
accuracy than the invention, but this is at the expense of using one receiver per antenna 
instead of one receiver only for all antennas. 



